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A precursor effect on the Fermi surface in the two-dimensional Hubbard model at finite tempera- 
tures near the antiferromagnetic instability is studied using three different itinerant approaches: the 
second order perturbation theory, the paramagnon theory (PT), and the two-particle self-consistent 
(TPSC) approach. In general, at finite temperature, the Fermi surface of the interacting electron 
systems is not sharply defined due to the broadening effects of the self-energy. In order to take ac- 
count of those effects we consider the single-particle spectral function v4(k, 0) at the Fermi level, to 
describe the counterpart of the Fermi surface at T = 0. We find that the Fermi surface is destroyed 
close to the pseudogap regime due to the spin-fluctuation effects in both PT and TPSC approaches. 
Moreover, the top of the effective valence band is located around k = (tt/2, 7t/2) in agreement with 
earlier investigations on the single-hole motion in the antiferromagnetic background. A crossover 
behavior from the Fermi-liquid regime to the pseudogap regime is observed in the electron concen- 
tration dependence of the spectral function and the self-energy. 

PACS numbers: 



I. INTRODUCTION 

Recently, investigations of the Fermi surface (FS) have been a crucial issue in the study on the high-T c cuprate 
superconductors. Angle- resolved photoemission experiments |l],^| have found evidences of deviations from conventional 
Fermi liquid demonstrating the existence of the pseudogap formation above T c in the single-particle spectral function 
and the breakdown of the FS at the so-called "hot-spot" region in the Brillouin zone. 

The Fermi surface is an important physical property which characterizes the low temperature behavior of the normal 
phase of fermion systems. If the stability of the Fermi liquid state is verified, the Luttinger theorem || concerning 
the conservation of the Fermi surface volume holds. However, if the interacting fermion system is in a state close to 
some instability, the precise definition of both the Fermi liquid and the FS may become uncertain. 

Some indications of the deviations from the conventional Fermi surface behavior have already been seen in the 
two-dimensional (2D) Hubbard model. Quantum Monte Carlo (QMC) simulations j| |) indicated the change of the 
FS topology for small hole dopings in the strong coupling regime. This topological change is characterized by the FS 
closing around k = (tt,tt) despite the (small) hole dopings. In the weak coupling regime of the Hubbard model, using 
second order perturbation theory, Zlatic, Schotte, and Schliecker obtained that the FS changes its topology due 
to the self-energy effects for small dopings. In contradiction with that, within the same framework of perturbation 
theory, Halboth and Metzner || have shown that the modification of the FS is very small for weak coupling at zero 
temperature. Thus, even in the weak-coupling regime, the FS of the Hubbard model leads to controversy. 

In this paper, we focus on the modification suffered by the FS of the 2D Hubbard model near the antiferromagnetic 
instability in the weak-coupling regime. The 2D Hubbard model produces several different electronic phases according 
to the choice of the parameters. For the highly doped case, the system is in the conventional Fermi liquid regime with 
a well-defined FS. On the other hand, in the ground state at half-filling, the system becomes an antiferromagnetic 
insulator. Thus, for an appropriate set of parameters, the zero temperature phase transition will take place, and quite 
probably, the conventional Fermi surface picture is lost before that. 
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However the long range antiferromagnetic order at half-filling is automatically destroyed at sufficiently high temper- 
atures or finite doping. Thus, assuming paramagnetic background at finite temperature we use itinerant approaches 
such as the second order perturbation theory, the paramagnon theory || and also the recently developed two-particle 
self-consistent approach JTo} |T^| to calculate the effects produced on the FS. Varying the electron concentration, the 
coupling constant and the temperature, the single-particle spectral function is calculated in all those approaches. In 
general the FS is well defined for a normal metal at T = 0. At finite temperature, however, we cannot define the FS 
sharply. To describe the counterpart of the Fermi surface at finite temperatures, we introduce the spectral function 
A(k, 0) at the Fermi level. The meaning of this quantity is discussed in detail. 

We will see that the FS of the Hubbard model at finite temperature shows remarkable modifications. As we 
approach the antiferromagnetic instability condition, the FS is destroyed at the (tt, 0) regions in the Brillouin zone. 
The calculated effective band dispersion indicates that the anisotropic pseudogap and the destruction of the FS are 
both manifestations of the precursor of the antiferromagnetic instability at half-filling. Besides that, we note that 
there is a trace of the FS in the pseudogap regime in which there is no effective band dispersion crossing the Fermi 
level. Also, by varying the electron number per site n, we find that the spectral function and self-energy clearly show 
a crossover from the Fermi liquid regime at the smaller n to the pseudogap regime near the half-filling. 



II. MODEL 



We start with the standard 2D Hubbard Hamiltonian defined by 

H = -t 22 a\ a a ja + U ^riijni^- ii 2_,n ifa (1) 

<i,j>,(T i i.cr 

where t is the nearest-neighbor transfer matrix on the square lattice, U is the on-site repulsive coupling constant, fj, 
the chemical potential and m a — al a ai a . The Fermi surface can be determined using the single-particle spectra which 
is calculated with the full Green's function. The spectral function A(k, u>) is given in terms of the imaginary part of 
the retarded Green's function: A(k, ui) = — (l/7r)IruG(k, ui) where G(k, ui) is given by 

G(k,tj)- 1 =w + 2 7 ? -(e k - M )-£(k,w) (2) 

with £k = — 2t(cosk x + cosfcj,) and £(k,a->) is a proper self-energy which is given in the following subsection. The 
chemical potential fj, is determined to be consistent with a given electron concentration n which satisfies 

/oo 
dujf(uj)A(k,uj) (3) 

for fixed values of U/t and T/t. Here /(u>) is the Fermi distribution function defined by f(co) = l/[exp(ui/T) + 1]. 

Next, we introduce self-energies for three different approaches; the paramagnon theory |9|], the second order per- 
turbation theory, and the two-particle self-consistent approach |Tc| . 



A. Paramagnon theory and second order perturbation theory 

The paramagnon theory (PT) is probably the simplest way to take account of the magnetic instability in the 
electronic states qualitatively. In that approach, the divergence of the spin susceptibility in the random phase ap- 
proximation (RPA) signals the antiferromagnetic instability. At T = 0, due to the perfect nesting condition, the real 
part of the one-loop polarization bubble RexoCcijO) diverges, and the PT approach cannot be simply applied, xa is 
given by 

Xo(,q, v m ) - > ; -. r (4J 

with v m — 2nnrT and fj, is the chemical potential in the non-interacting case. At non-zero temperature, however, 
Re%o(q, 0) remains finite even at half- filling. Thus, the RPA spin susceptibility does not diverge for suitable small 
values of the coupling constant U and the system remains in its paramagnetic regime. We perform all calculations in 
the present paper with suitable choices of the temperature T, the coupling constant U, and the electron concentration 
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n. A zero-temperature version of our approximation has been used by Kampf to calculate the spectral function 
of the doped Hubbard model. One drawback of the paramagnon theory is that it does predict the antiferromagnetic 
instability in two dimensions violating in this way the Mermin- Wagner theorem. However, if we avoid this transi- 
tion regime by working at finite temperatures, as we have shown in our previous work ||, this approach gives the 
qualitatively correct non-Fermi-liquid energy dependence in the self-energy with the corresponding anisotropic pseu- 
dogap in the single-particle spectral function. These results are in general agreement with several other approaches 



P-l| , |l2 14 ljj. Moreover, the anisotropy of the pseudogap, which has been observed also in the QMC simulation by 
Crcffield et al |l(|, was shown in our previous paper. 

To derive the self-energy, we use the standard diagram technique by treating the interaction term perturbatively. 
Since the Hartree-Fock term Un/2 is nothing but a constant energy shift, we include it in the chemical potential. The 
self-energy in the PT framework can be written as 



E PT (k, W „) = [/ 2 T]T£ 
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-Xo(q, fm)]Gti(k + q,w n + z/ m ) (5) 

where C?o(k, w n ) is the non-interacting Green's function with u> n = (2n + 1)tvT. Here we have introduced the charge 
susceptibility by Xcfavm) = Xo(q, v m )/[l + Uxo(<l,Vm)] and the spin susceptibility by Xs(l,^m) = Xo(q,^m)/[l - 
Uxa(<l, I'm)]- In the angular bracket in Eq. (||), Xoi^^m) is subtracted to avoid the double counting of the lowest 
order diagram. £(k, uS) in Eq. (||) is obtained from E PT (k, ui n ) by performing an analytic continuation. 
In contrast, the self-energy obtained from conventional second order perturbation theory is given by 

S 2nd (k, u n ) = U 2 TY,Y1 *o(q, OG (k + q, u n + v m ). (6) 

Urn q 

This is equivalent to the lowest order contribution of the self-energy S PT (k,o;, i ) in the paramagnon theory. Detailed 
investigation on S 2nd (k, u> n ) of the Hubbard model is seen in references p|jqjr^|. 



B. Two-particle self-consistent approach 

The two-particle self-consistent (TPSC) approach has been developed by Vilk and Tremblay et al [p~0| p~2|| - This 
approach can give both the spin and charge structure factors in good quantitative agreement with the results of the 
quantum Monte Carlo simulations in the Hubbard model Jlcfl . Detailed argument on this approach has been given in 
their review article [^2| . Following their formulation, the self-energy of the single electron Green's function in TPSC 
approach is given by 

E TPSC (k,^) = ^tJ2J2 [twr sc (q,M 

v m q 

+£4 P xJ PSC (q, v m )] G (k + q, u n + u m ) (7) 

where X<f PSC (q, v m ) = xo(q, v m )/[l + C/ c hXo(q, v m )\ and xJ PS °(q, v m ) = Xo(q, v m )/[l - C/ sp Xo(q, v m )\. The Hartree- 
Fock term has been absorbed again in the chemical potential. Note the difference of the factor 2 of our definition of 
Xo here from the one in their original papers. The renormalized coupling constants U sp and t/ c h are determined by 
solving the self-consistent equations given by 

rEE2 X J PSC (q,, m ) = n -2( Ilt r il ), (8) 

Km q 



T E E 2 Xc raU (q, v m) = n + 2 (n T n x ) - n\ (9) 

K m q 

and 

4C7 

U ** = ^2 ^ n i) ■ ( 10 ) 
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Here (n-frij.) is the average double occupancy per site. 

Note that contrary to the other two methods the Mermin- Wagner theorem is satisfied at any finite temperature in 
the TPSC approach To review it briefly let's consider the half-filling case. At half-filling, from Eqs. and 

(^), we obtain the following self-consistent equation to compute U sp : 

T W 2xo(q,^ m ) = U sp . , 

If U sp always satisfies the condition 1 > C/ sp xo(Q,0), the magnetic instability is not realized, namley the Mermin- 
Wagner theorem holds. It is convenient to rewrite the above self-consistent equation in terms of two functions fi(u) 
and fziu) defined by 



/i(ti) = rVV 2 *°(y™) ( i2) 



and 



h(u) = l-^j. (13) 

The self-consistent equation ( [TT| ) is equivalent to the relation fi(U sp ) = f2(U sp ). Fig. [j] shows a schematic plot of 
fi(u) and f2(u) as the functions of u for finite T. The behavior of f2{u) is trivial. Note that the slope of /^(n) 
is in proportion with X/U. Using the sum rule of xo given in Appendix, we see fi(u = 0) = 0.5. diverges 
as u approaches l/xo(Q,0). This divergence was proved by Vilk and Tremblay |l2). The intersection of fi(u) and 
/2(it) gives the solution of the self-consistent equation. From the plot, we see that the intersection point always stays 
within the range < u < l/xo(Q> 0) for arbitrary U. Namely, the solution of the TPSC self-consistent equation never 
violates the Mermin- Wagner theorem. 



III. FERMI SURFACE AT FINITE TEMPERATURE 

Let's consider how to obtain the information on the Fermi surface from the full Green's function given by Eq. (||) 
at finite temperature. So far, several different methods to investigate the FS have been used in the literature. One 
straightforward way to define the FS is given by the quasiparticle poles of the Green's function through the solutions 
of £k F — A* + ReE(ki?,0) = 0. Another way is to use the electron momentum distribution function u(k) defined by 
n(k) = f_ du>f(uj)A(k, uj). In this case, the FS is characterized by each k point associated with a finite discontinuity 
of n(k). Those definitions give a proper FS at zero temperature. However, at finite temperature, since the imaginary 
part of the sclf-cncrgy has a finite value with a non-trivial k dependence in general, these definitions are not applicable. 
Another definition is to use n(kp) — 1/2. This is not useful to apply for the pseudogap regime since this condition is 
always satisfied at half-filling and at finite temperature. Maximum slope of Vk«(k) is also used in the determination 
of the FS obtained with the angle-resolved photoemission |l8|,[l9|]. However, because 77,(k) is obtained from a quantity 
which is integrated over the energy, it has a maximum slope even if a gap or a pseudogap exists. At finite T the FS 
is often defined approximately by the k points at which the effective band dispersion crosses on the Fermi level. The 
effective band dispersion is obtained from the maximum intensity points of the spectral function. 

Strictly speaking, at finite temperature we cannot define a sharp FS because of the broadening effects in the 
spectral function. However, in principle, the Fermi level can be calculated for any sets of parameters of the system if 
the number of electrons is known. We therefore use this route to establish the Fermi surface. Taking into account the 
above consideration on the FS, instead of employing conventional definitions, we consider the electron distribution at 
the Fermi level. Such distribution can be represented by a single-particle spectral function written as 

ACk 0) = 1 |ImS(k,0)| (u 

V ' ' ^[£ k -/!-Rc£(k,0)] 2 + [ImE(k,0)] 2 ' 1 ' 

It is evident that it satisfies a sum-rule N(Q) = ^ k A(k, 0) where N(0) is the density of states (DOS) at the Fermi level. 
A(k, 0) contains all the information of the single-particle states at the Fermi level and possesses enough information 
to study the properties of the Fermi surface. 

We easily see that A(k, 0) has the following features, which covers the basic properties of the conventional FS. 
When the system is described by the conventional Fermi-liquid at zero temperature, we have that 
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A(k, 0) = S(e k - [i- Re£(k, 0)) 



(15) 



and the peak points of the delta-function trace a sharply defined FS in the k space. Thus, A(k, 0) recovers the 
conventional FS at the ground states. In a Fermi liquid at finite temperature, the spectral function becomes broad 
around the Fermi level. Quite naturally, A(k, 0) covers this feature because it has a finite width and it draws a broad 
trace in the k space. On the other hand, when the system does not have a spectral weight at all at the Fermi level 
(N(Q) — 0) - for instance, an insulating state with a real gap - A(k, 0) is identically zero for all wave vectors. Thus, 
A(k, 0) represents all cases ranging from the well-defined FS to the vanishing FS limit. As a result, A(k, 0) can be 
used to trace the evolution of the FS in all cases. Observed Fermi surface in angle-resolved photoemission experiments 
is extracted from the finite weight of the single particle spectra in the vicinity of the Fermi level. Here too, in this 
respect, it is meaningful to introduce A(k, 0) in order to study the FS structure. 

As we see from the argument above, the generalized FS at finite temperatures is no longer strictly speaking a surface 
in the momentum space. It is not clear if it could be used as the counterpart of the T = FS. This issue could be 
solved by analogy with what happens to the band dispersion in both non-interacting and interacting systems. The 
electronic states in the free-fermion system in a lattice is completely described by a non-interacting band dispersion 
lu = E k — [i. In the interacting system, within certain theoretical treatment of the interaction term, the spectral 
function A(k, to) contains all the information of the single-particle states, but the (effective) band dispersion cannot 
fully describe the electronic states. The FS of the non-interacting system is given by E k — /x = 0. On the other hand, 
A(k, 0) at the Fermi level continues to carry all necessary information concerning the electronic states. As a result 
it replaces the band dispersion in the description of the electronic states in the interacting system. In the following 
section of this paper, we use the term, Fermi surface, in this extended sense. 



IV. RESULTS 



Using the self-energies and the Green's functions presented in the previous section, numerical calculations were 
performed in the same manner as in our previous work || . A mesh of about 4000 grids in the first Brillouin zone has 
been used to compute the momentum integral to obtain the imaginary parts of both charge and spin susceptibilities. 
A 120 x 120 mesh on the Brillouin zone has been used to perform the momentum integrations for the self-energies to 
obtain the Green's functions and the spectral functions. 



A. Second order perturbation theory 

In Fig. |^ we show the temperature dependence of DOS calculated using the second order perturbation theory for 
U = 4.0i and n — 1.0. Detailed analyses of the self-energy and the single-particle spectra at finite temperatures 
and finite dopings have already been done by Zlatic, Grabowski, and Entel |17|] . A main temperature dependence 
in DOS appears only around the Fermi level (u = 0). We do not see any pseudogap-like structure in DOS in this 
approximation. At T = 0.5i, the maximum of DOS at the Fermi level is broad. For a lower temperature T = O.lt, 
this maximum becomes a sharp peak. 

Figure || shows the gray-scale density plot of A(k, 0) for the corresponding parameters used in Fig. 0. As emphasized 
earlier on, the density plot of A(k, 0) in Fig. || gives a counterpart of the Fermi surface. In this plot, the larger (smaller 
or zero) intensity of A(k, 0) is represented by the bright (dark) region and the axes for both k x and k y directions are 
normalized by 7T. For T = O.lt (lower panel), a well-defined FS has been obtained. The broadening in the plot for 
T = Q.ht (upper panel) can be understood as the temperature effect in the self-energy. Thus, within our results, the 
second order perturbation theory gives the conventional Fermi-liquid-like behavior both in the single-particle spectra 
and in the FS. Namely, this approach does not explain the breakdown of the Fermi surface. Our calculation has 



been done at half-filling only whereas the calculation of Zlatic et al |17| is for large doping (n = 0.8). They found a 



topology change of FS, but their choice of parameters may not be entirely suitable for that regime. 



B. Paramagnon theory 

Figure |^ shows the n-dependence of A(k, 0) for U — 2.0t at T = 0.22i calculated in terms of S PT (k, u>„). For 
n = 0.8 (upper panel), a clear and large FS is obtained indicating that the system is in the conventional Fermi liquid 
phase. For n — 0.9 (middle panel), the FS volume is enlarged by the increment of the electron concentration. In 
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addition, the distribution around (tt, 0) becomes broader. For n = 1.0 (lower panel), the FS around (ir,0) regions 
is completely destroyed and there remains well-defined Fermi surface segments only around (7r/2, 7r/2) regions. This 
anisotropic destruction of the FS is caused by the anisotropy of the pseudogap formation. The origin of this kind 
of anisotropic pseudogap has been explained in our previous work The strong low-energy antiferromagnetic 
Stoner enhancement in the spin susceptibility produces the non-Fermi-liquid energy dependence in the fermion self- 
energy. The momentum dependence of the self-energy is determined by the anisotropy of the band flatness of the 
non-interacting band dispersion at the Fermi level. As a result of those effects, the pseudogap opens first at the (tt, 0) 
regions, and the FS is partially destroyed. 

In Fig. |B| we plot the [/-dependence of the spectral function A(k, 0) within the paramagnon theory for T = 0.5i 
at half-filling. For U = 2.0t (upper panel), the FS is rather sharply defined as for the non-interacting case. As U 
becomes large, the distribution becomes broad and the FS cannot be sharply defined. The tendency for broadening 
is remarkable in the (tt, 0) regions. The original structure of the spectral function around (tt, 0) is almost washed out 
for U = 3.0t. Instead, there is a hollow structure centered on (ir,0). The (tt/2,tt/2) regions have relatively large 
intensities. 

In Fig. ^| we show the effective band dispersion obtained from the main peaks of the spectral functions A(k, to) for 
the same parameter as the lower panel in Fig. ^. The pseudogap around (tt, 0) is larger than the one at (tt/2, 7r/2), and 
segments of the shadow bands ]|(| can be seen around this point. This result indicates that the top of the effective 
lower (valence) band is located at k = (ir/2,ir/2). This result is consistent with the earlier investigations of the 
single- hole motion in the 2D antiferromagnetic background f2~]]j . In Fig. ^ it is important to realize that the effective 
band does not cross the Fermi level. In conventional meaning we cannot define the FS in this regime, but the spectral 
weight remains at the Fermi level and it forms a very broad remnant Fermi surface structure in the momentum space 
as we see in the lower-most panel of Fig. [5l 



C. Two-particle self-consistent approach 

Figure | shows the plot of A(k, 0) in the TPSC for U = 4.0t at half-filling for T = (upper panel) and T = QM 
(lower panel). Despite the above mentioned difference between the TPSC and PT, both approaches show qualitetive 
agreement with respect to the destruction of FS. The quantitative diffference between the two methods is due to the 
fact that in the PT case the interaction effects are overestimated by the lack of vertex corrections. 

In Fig. U we show the temperature dependence of DOS N(to) in the TPSC approach for U — At and n = 1. As we 
lower the temperature, the gap-like structure grows at the Fermi energy. That behavior is well correlated with the 
destruction of the FS observed in Fig. ^. 

The electron concentration dependence of the spectral intensity A(k, 0) is shown in Fig. [)] for U — At and T = 0.25i. 
For n = 0.85 (upper panel), the FS is still defined well. At n — 0.9, the FS enlarged but in the regions near the 
Brillouin zone it gets broader. At n = 0.95, the FS region near (tt, 0) becomes very unclear, while the other FS sectors 
remain relatively well-defined. Again, this indicates the anisotropic breakdown of the Fermi liquid picture on the 
Fermi level. 



In Fig. 10 we plot the spectral functions and the self-energies for various electron concentrations for T = 0.25i and 
U = AM. For convenience, the momentum kp has been assumed to be the same one in the (1,0) direction on the 
non-interacting Fermi surfcae for each n. At n — 0.8, the spectral function (a) has a single peak at the Fermi level. 
This peak corresponds to the conventional Fermi liquid quasiparticle. As increasing n towards half-filling, the peak 
becomes broader and shifts into the u> < side. On the other hand, a shoulder appears in the u> > side in the 
spectral function (n — 0.9) and in the end this turns into another peak. Thus the quasiparticle is destroyed, and 
two-peak structure emerges at n = 1. The self-energies also show the breakdown of the Fermi liquid picture. In our 
standard knowledge on the Fermi liquid theory, the real part of the self-energy ReE is linear in to and the slope is 
negative, and the imaginary part of the self-energy ImS « w 2 . At n = 0.8, the real (b) and imaginary (c) parts of 
the self-energy show those Fermi liquid like behaviors. As n becomes larger, the self-energy dramatically changes its 
lo dependence and indicates deviations from the Fermi liquid nature. There are two remarkable features at n = 1: (i) 
the lo dependence of the real part of the self-energy resembles to the spin-density- wave self-energy (~ 1 /to), (ii) the 
large negative peak appears in the imaginary part of the self-energy. 
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V. DISCUSSION AND CONCLUSION 



We find that the second order perturbation theory fails to show the breakdown of the Fermi surface. The decisive 
difference between the second order perturbation and the other two approaches is the fact that in the former case the 
Stoner factor is not taken into account in the self-energy. We can immediately see this by comparing the definitions 
of these self-energies. As we have already mentioned, the Stoner enhancement is important to produce the anomalous 
energy dependence in the self-energy. If it is combined with the flat-band dispersion the anisotropic pseudogap is 
formed in the spectral function. Because of this difference in self-energy, the second order perturbation theory FS is 
well defined and does not show any anomaly even for dopings close to half-filling. 

As mentioned in Section 1 there has been some controversy concerning the Fermi surface topology in the 2nd 
order perturbation theory framework j?],||,[l7]]. This issue was recently discussed by Halboth and Metzner (HM) in 
||. We just want to add that in our calculation we pay attention to the existence of different values of chemical 
potentials for the non-interacting Green's function and the interacting Green's functions which is renormalized by the 
self-energy effects. If one uses a fixed value of chemical potential in both cases the conservation of electron number 
is automatically destroyed. This will in turn produce discrepancies for any finite doping. This does not seem to be 
taken into account in Refs. R] and Jl7| . 

We discuss next the possibility of the small hole pocket formation in the 2D Hubbard model. Langer et al 0] and 
Schmalianef al [23f| , using the fluctuation exchange (FLEX) approximation for the Hubbard model, have investigated 
how the low-energy electronic states and the FS evolve at finite dopings. They explained |Q one scenario of the 
formation of the small hole pockets produced by the transfer of the spectral weight from the main band to the shadow 
band. In the range of the parameters used in our calculation, we have not obtained any sign of the hole pockets or 
even the "shadow" of the Fermi surface 0| , while we observed the segments of the shadow bands in the effective 
band dispersion. According to the argument in Ref. f22| , the origin of the shadow band is a large imaginary part of 
the self-energy at the k + Q regions rather than another pole in the single-particle Green's function. In their result 
the imaginary part of the self-energy has a strong momentum dependence. That feature has been seen also in our 
previous calculation ||. However, the momentum dependence of ImE(k, lo) produces the anisotropy of the pseudogap 
in A(k, oj). In both PT and TPSC approach, the origin of the double-peak is clearly associated with the real poles of 
the full Green's function. Because of this difference, it is natural that our data showed neither the shadow of the FS 
nor the hole-pockets. 

Chubukov and Morr 0] have shown the small hole pocket formation in an analytical investigation based on the 
nearly antiferromagnetic Fermi liquid approach for the spin-fermion model. Schmalian et al |l5|] have also obtained the 
small half-pocket-like segments in the FS for the spin-fermion model. Although both the PT and TPSC approaches 
are in a different approximation level from their semi-phenomenological model, the importance of the spin-fluctuation 
is common in both approaches. To compare with their results we need to study the t-t'-U Hubbard model p5| which 
has the same non-interacting band dispersion including the next-nearest-neighbor hopping t' and to use the Fermi 
surface topology as they used. If the same mechanism discussed in the present paper works in the t-t'-U model, we 
can expect to obtain a half closed small FS around the (tt/2,it/2) point. 

The single-particle spectrum of the t-J model is equivalent to that of the Hubbard model in the strong coupling 
limit if we take no consideration of the upper Hubbard band effect. Recently, Dai and Su |^6| obtained the destruction 
of the FS of the t-J model, and the FS for a hole doping clearly shows a topological change of the Fermi surface which 
is consistent with the results of the QMC simulations in the large-coupling Hubbard model jl| ||. Putikka, Luchini, 
and Singh [^7j also showed the violation of the Luttinger theorem in the t-J model. 

In the present work which deals with the weak-coupling Hubbard model, we cannot make any decisive statement 
concerning the validity of the Luttinger theorem although the complete FS structure no longer exists. At half- filling 
in the arbitrary coupling, the Hubbard model acquires the particle-hole symmetry. Hence, presumably the half-filling 
Fermi surface (if it exists) cannot show any anomalous topology. As it is well known, in the large coupling regime, 
the doped Hubbard model have a preformed gap associated with the Mott-Hubbard transition. The energy scale 
of the preformed gap is of the order of U. The QMC simulations have shown that the doped hole forms a narrow 
quasiparticle band dispersion on the top of the lower Hubbard band Namely, there is a jump of the Fermi level 

between the infinitesimal hole and electron dopings. This discontinuity of the Fermi level produces the deviation of 
the electron-hole symmetry in the electronic states. In the weak-coupling regime, since the effect of the preformed 
gap at finite temperatures is considerably weak, the electronic states must preserve the electron-hole symmetry for 
an infinitesimal doping, and the topological change cannot be seen. Thus, the topology change of the FS in the 
large coupling regime may be a manifestation of the effective violation of the electron-hole symmetry due to the large 
preformed gap. 
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In conclusion, we studied the electron distribution A(k, 0) at the Fermi level for the 2D Hubbard model in the 
weak-coupling regime. This quantity is useful for the generalization of the conventional Fermi surface concept and 
can be applied in all circumstances for all ranges of parameters of the Hubbard model at finite temperature. We 
examined three different itinerant approaches, and found that the second order perturbation theory fails to describe 
the breakdown of the Fermi surface. In the pseudogap regime, we found that the FS is partially destroyed in PT 
and TPSC. From the anisotropic pseudogap and also from the destruction of the FS, we observe that the top of the 
effective valence band is located at k = (tt/2,tt/2). This is in general agreement with prevoius theoretical studies for 
a strongly coupled regime |2l]] . We also investigated the u dependence of the spectral function and the self-energy on 
the FS. We showed that there exists only quasiparticle like state for small n. In contrast, as we increase n, we observe 
a crossover between a Fermi-liquid phase and a spin-density-wave-like regime with the evolution of the precursor 
pseudogap. 
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APPENDIX: A SUM RULE FOR xo(Q,^m) 

In this appendix, we introduce a sum rule for xo(<L v m) given by 

n /_ n 



r££xo(q,*») = 5(i-5). (Ai) 

v m q 

This holds for arbitrary electron concentrations and for any temperature. It can be proved by calling the definition 
of xo(qj v m ) given by Eq. (Q). Physically speaking, this sum rule states that the total summation of the particle- hole 
bubble Xo(q> v m ) over q (in the Brillouin zone) and v m is equal to the products of the electron concentration n/2 and 
the hole concentration 1 — {n/2) per spin. It is possible to generalize this sum rule for the one-loop diagram given in 
terms of the full Green's functions G(k, u> n ). In this case, the sum rule becomes 

r££x(q,^) = |(l-|) (A2) 
v m q 

where 



It can be proved in terms of 



X(q,^m) = -T££G(k + q,w„ + v m )G(k,Lu n ). (A3) 

u, k 



, duj' ImG(k,a/) . . 

G(k,w„) = / , . - ■ (A4) 

, 7T U! — lUJ n 
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FIG. 1. A schematic plot of the functions fi(u) and /2(H). The solid circles A and B indicate the solutions of the equation 
fi(u) = /2(u) for larger U and smaller U, respectively. 

FIG. 2. Temperature dependence of the density of states N(co) within the second order perturbation theory for U = 4.04 at 
half-filling (n = 1). 



FIG. 3. Spectral intensity A(k, 0) in the Brillouin zone calculated within the second order perturbation for U = 4.04 at 
half-filling (n = 1). Upper panel for T = 0.54 and lower panel for T = 0.14. 



FIG. 4. Spectral intensity A(k, 0) in the Brillouin zone within the paramagnon theory calculated for various electron 
concentration n for U = 2.04 at T = 0.224. The upper panel is for n = 0.8, the middle panel for n = 0.9, and the bottom panel 
for n = 1.0. 



FIG. 5. Spectral intensity A(k, 0) in the Brillouin zone within the paramagnon theory calculated for various values of U for 
T = 0.5 at half-filling. From top to bottom, U = 2.04, 2.54 and 3.04, respectively. 



FIG. 6. Effective band dispersion obtained from the spectral function A(k,uj) within the paramagnon theory calculated for 
U = 3.04, T = 0.54 at half-filling. 

FIG. 7. Temperature dependence of the spectral intensity A(k, 0) in the Brillouin zone calculated within the two-particle 
self-consistent approach for U = 4.04 and n = 1. Upper panel for T = 0.54 and lower panel for T = 0.34. 
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FIG. 8. Temperature dependence of the density of states calculated within the two-particle self-consistent approach for 
T — 0.3t (solid curve) and 0.5t (broken curve) at U = 4t and n = 1. 

FIG. 9. Electron concentration n dependence of the spectral intensity A(k, 0) in the Brillouin zone calculated within the 
two-particle self-consistent approach for U = 4.0t and T = 0.25. The upper panel is for n — 0.85, the middle panel for n = 0.9, 
and the lower panel for n = 0.95. 

FIG. 10. The spectral function (a), the real part (b) and imaginary part (c) of the self-energy within the two-particle 
self-consistent approach at U = At and T = 0.25t for various electron concentration n. The Fermi momentum in the 
(1, 0)-direction on the non-interacting Fermi surface has been used as for each n. 
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